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ABSTRACT 

The radio source Hercules A resides at the center of a cooling flow cluster of galaxies at redshift z = 0.154. 
A Chandra X-ray image reveals a shock front in the intracluster medium (ICM) surrounding the radio source, 
about 160 kpc from the active galactic nucleus (AGN) that hosts it. The shock has a Mach number of 1.65, 
making it the strongest of the cluster-scale shocks driven by an AGN outburst found so far. The age of the 
outburst ~ 5.9 x 10 7 y, its energy ~3x 10 61 erg and its mean power ~ 1.6 x 10 46 erg s" 1 . As for the other 
large AGN outbursts in cooling flow clusters, this outburst overwhelms radiative losses from the ICM of the 
Hercules A cluster by a factor of ~ 100. It adds to the case that AGN outbursts are a significant source of 
preheating for the ICM. Unless the mechanical efficiency of the AGN in Hercules A exceeds 10%, the central 
black hole must have grown by more than 1 .7 x 10 8 M Q to power this one outburst. 

Subject headings: cooling flows - galaxies: clusters: individual (Hercules A) - intergalactic medium - X-rays: 
galaxies: clusters 



1. INTRODUCTION 

It is not yet clear which heating mec hanism (e.g. 
iNaravan & Medvedev"ll2001t iMotl et al.ll2004 is chiefly re- 
sponsi ble for preventing ga s from cooling in cluster cooling 
flows (Peterson et al. 2003), b ut the most promising is heat- 
ing by a central AGN JTabor & Binnevll993llTucker ■ & Davidl 
1997). Heating and cooling rates are linked if the AGN is fed 
by cooled or cooling gas. Such feedback could maintain oth- 
erwise unstable co ol cor es, explaining the prevalence of cool- 
ing flows (lFa bian|ll9 94T). AGN powered radio lobe cavities 
(e.g. Carill i et aljll994t IMcNamara et alJl2000t iFabian et al.1 
2000) heat the ICM JChurazoy et al.l2002l) but n ot enough to 
make up for radiative losses (Birzan et al. 2004). Weak "co- 
coon " shocks, long expected in model s of jet-fed radio lobes 
(e.g. IScheuerlll974l iHeinz et al]ll?98l). have been found in a 
numb er of system s (Fabian et al. 2003: iForman et aTl 120051 
IMcNamara et al.l 120051 iNulsen et al.l 120051) . They repre sent 
additional heating due to AGN outbursts, and provide a new 
tool for determining ages and energies of AGN outbursts. 

This letter reports the discovery of a shock front gener- 
ated by the AGN outburst that powers H ercules A. One of the 
bright est radio sources in the sky JDreher & Feigelsoril ll984l 
iGizani &L eahv 2003), Hercules A resides at the center of a 
cluster of galaxies with X-ray luminosity ~ 5 x 10 44 erg s" 1 , 
at red shift z = 0.154 iSiebertet al.l Il999t IGizani & Leahvl 
2004). Despite its high radio luminosity, Hercules A lacks 
bright radio hotspots and so belongs to Fanaroff-Riley class 
I, but with an unusual, j et-dominate d, radio m orphology 
JDreher & Feigelsonl Il984t IGizani & Leahvl I20031) . Using 
Einstein spectra, iWhite et al.1 lJl99l) found a formal cooling 
rate of zero for the Hercules A Cluster, but the high central 
density we find (n e > 0.1 cm" 3 ) gives it a central cooling time 
typical of a cooling flow cluster. 

Section |2]gives details of the observations and data reduc- 
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tion, and section|3]discusses the main features of the Chandra 
image of Hercules A. Properties of the shock are discussed in 
section|4]and its implications in section|5] Flat ACDM, with 
Hq = 70 km s" 1 Mpc -1 and 57 m = 0.3, is assumed throughout, 
giving a scale of 2.67 kpc arcsec" 1 for Hercules A. 

2. OBSERVATIONS AND DATA REDUCTION 

Hercules A was observed with Chandra for 14.8 ksec on 
25 Jul 2001, in VFAINT mode with ACIS-S at the aim point 
(OBSID 1625). For the analysis here, the event list was repro- 
cessed using recent calibrations. It was screened to remove 
ASCA grades 1, 5 and 7, and bad pixels. Periods of high 
particle background were removed following the method of 
Markevitch 6 , leaving 12.4 ksec of good exposure time. After 
cleaning, the mean count rate in ACIS SI was 0.136 ct s , 
~ 3.5er (9%) higher than expected 7 , suggesting some residual 
contamination due to particle background. Data were pro- 
cessed to correct for time dependence of the ACIS gain 8 and 
filtered according to the prescription of Vikhlinin 9 to reduce 
particle background. Background event files were created by 
processing standard ACIS background files in the same man- 
ner as the data. Point sources were identified manually for 
removal from spectra and surface brightness profiles. ARF's 
and RMF's for extended regions are weighted by number of 
events. ARF's are corrected to allow for the reduction in low 
energy response due to contaminant on the ACIS filters. 

3. THE X-RAY IMAGE OF HERCULES A 

The image in the upper panel of Fig. ^ shows the filtered 
and calibrated events from the raw Chandra data (evt2 file) for 
Hercules A, binned by a factor of 4. The lower panel shows 
a 0.3 - 7.5 keV image made after the cleaning and reprocess- 
ing described above. The image has been smoothed with a 2" 
gaussian and divided by a beta model, with 42" core radius 
and a beta of 0.6. cen tered on the X-ray peak (core radius from 
Gizani & Leahv l2"o"04l but smaller beta). Division by the beta 
model reduces the radial variation of surface brightness, mak- 

6 http://hea-www.harvard.edu/ maxim/axaf/acisbg 

7 http://hea-www.harvard.edu/ maxim/axaf/acisbg/data/README 

8 http://hea-www.harvard.edu/~ alexey/acis/tgain/ 

9 http://cxc.harvard.edu/cal/Acis/Cal_prods/vfbkgrnd/ 
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FIG. 1 . — Upper panel: 0.3 - 7.5 keV Chandra image of Hercules A made 
from the distributed evt2 file binned by a factor of 4. Lower panel: 0.3 - 7.5 
keV image of Hercules A made from the cleaned, reprocessed data, smoothed 
with a 2" gaussian and divided by a beta model. The scale bar in each panel is 
1 arcmin (160 kpc) in length. The bright central region ~ 1 arcmin in radius 
is surrounded by the shock front. The southwest cavity is ~ 0'.5 from the 
central peak. 



ing it easier to discern substructure over a substantial range of 
radius. Each image has a 1' scale bar. 

Although the central peak of the X-ray image is prominent, 
it is well resol ved by Chandra and there is no sign of a point- 
like AGN Ccf. lTrussoni et alJl2001l) . A striking feature of the 
X-ray image is the bright region, roughly 1' in radius, that 
stands out in the upper panel of Fig.^ This has a similar size 
to the radio emission and extends to the east and west around 
the radio lobes (Fig. flower, Fig.|2ji. Its shape and association 
with the radio source suggest it is the shocked cocoon of the 
expanding radio lobes ( Scheuer 1974t lHeinz et alJl998l) . The 
break in surface brightness that bounds this region is shown 
to be consistent with a shock front below. 

There is a ~ la deficit of X-ray emission in the region ~ 
0'.5 to the southwest of the bright center in Fig. [2 ~ 15" (40 
kpc) in radius. There is a weaker, ~ 3a, deficit in the X- 
ray emission from the corresponding region to the northeast, 




FIG. 2. — X-ray and radio images of Hercules A. The Chandra imag e of 
Fig.Qoverlaid with 1.4 GHz radio contours from lGizam &Leah\ll2Tjrj3l) . 



partly masked by a bright spot of X-ray emission to the north 
of the center. These features resemble the c avities associated 
with many other cluster radio sources (e.g. iMcNamara et alJ 
2000; Fabi an et al.l20 00). However, they are not aligned with 
the axis of the radio jets an d do not contain radio l obes. They 
might be ghost cavities (e.g. McNamara et al. 2001), but if so, 
it is surprising that they lie within an active radio source. 

Lastly, there is a ridge of enhanced X-ray emission cross- 
ing the the bright region, from ~ 30° south of east to ~ 30° 
north of west, roughly at right angles to the axis defined by 
the cavities. This feature also has no obvious association with 
the radio source (it forms an angle of ~ 20° with the radio 
jets). The excess emission appears to be thermal, due to rel- 
atively cool, dense gas, which cannot be fully supported by 
hydrostatic forces. The gas may be cool filame nts, like those 
seen in other cluster cent ral radio sources (e.g. iForman et alJ 
120051 IrTulsen et ai1 l2005). or it may be a cooler disk that is 
partly supported by rotation. 

4. THE SHOCK FRONT IN HERCULES A 

The surface brightness profile of the bright circular region 
was measured in two 80° sectors, approximately at right an- 
gles to the axis of the radio jet. This avoids smearing the edge 
in the surface brightness profile due to elongation of the bright 
region to the east and west. Fig. [3] shows the radial surface 
brightness profile for the ranges of PA 330° - 50° and 150° - 
230° combined. Point sources were eliminated, background 
subtracted and the resulting profile exposure corrected. Al- 
though the data are quite noisy, there is a clear break in surface 
brightness at a radius of 60" (~ 160 kpc), at the edge of the 
bright central region. Beyond the break, the surface bright- 
ness is well fitted by the power law, r~ a with a = 2.1 1 ± 0.43 
(90%). We now consider the interpretation of this front as a 
shock. 

To determine the strength of the shock we use a spherically 
symmetric, hydrodynamic model of a point explosion at the 
center of an initially isothermal, hydrostatic atmosphere. Be- 
fore passage of the shock, the gas density is assumed to fol- 
low the power law, p(r) oc r"'', with r\ = 1 .55, chosen to make 
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FIG. 3. — Surface brightness profile of the shock front in Hercules A. The 
0.6 - 7.5 keV surface brightness profile is measured in sectors from PA 330° 
to 50° and 150° to 230°, to the north and south of the AGN, at right angles 
to the jet axis. Surface brightness errors are \a statistical errors. Radial error 
bars show the limits of the bins. The smooth curves are surface brightness 
profiles for shock models with Mach numbers of 1.51, 1.65 and 1.79, from 
bottom to top on the right. Models are scaled to match the observed surface 
brightness outside the shock. 



the surface brightness profile of the undisturbed gas match 
the observed profile outside the shock. The gravitational field 
(g oc 1 jr) is scaled to make the undisturbed atmosphere hy- 
drostatic. The surface brightness profile is determined from 
the model, assuming that the temperature of the unshocked 
gas is 4 keV (see below). Relative Chandra count rates in 
the 0.6 - 7.5 keV band are computed using detector response 
files from near to the aim point for these observations. The 
XSPEC wabs x mekal spectral model was used, with a fore- 
ground column density of 6.4 x 10 20 cm" 2 , a redshift of 0. 154 
and abundances of 0.5 times solar, appropriate for Hercules A 
(model surface brightness profiles are insensitive to these and 
the preshock temperature in the relevant temperature range). 
The model is self-similar, allowing it to be scaled in radius to 
match the location of the shock and in normalization to match 
observed surface brightness outside the shock. 

In Fig. [5] we show surface brightness profiles for model 
shocks with Mach numbers of 1.51, 1.65 and 1.79. A Mach 
1 .65 shock gives a reasonable fit to the data. Apart from the 
scaling, model parameters (the initial density power-law, r\, 
and preshock temperature) are constrained by observations, 
leaving only the Mach number of the shock free in the fit. 
The model has a number of shortcomings (the actual outburst 
is aspherical, does not inject energy in a single explosion and 
the initial gas density is not a power law, Nulsen et al. 2005), 
so that it can only be expected to match the data over a lim- 
ited range of radius behind the shock. Nevertheless, the fit 
provides a stringent test that this feature is due to a shock 
propagating into the cluster. 

In order to determine physical properties of the outburst 
from the model, we must determine the density and tempera- 
ture of the unshocked gas. However, outside the shock from 
1' to 2'.5, in the sectors of the surface brightness profile, there 
are only ~ 1250 photons in the 0.6 - 7.5 keV band. We have 
therefore used a single spectrum extracted from this region 
to determine the temperature and normalize the density of 
the unshocked gas. Using an absorbed mekal model, with 
Af H = 6.4 x 10 20 cm" 2 , redshift z = 0. 154 and the abundance set 
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FIG. 4. — Electron density profile of the shock front in Hercules A. Depro- 
jected electron density versus radius in the PA ranges 330° to 50° and 150° to 
230° . The shock is at 158 kpc. Density error bars are 90% confidence ranges. 
The continuous line shows the density profile for the best fitting model. 



to 0.5, gives a temperature of kT = 3.9^ ke V (90%). This is 
consi stent with previous measurements (e.g. Giz ani & Leahvl 
120031) . suggesting that the spectrum is not significantly af- 
fected by the particle background. Assuming that the gas is 
spherically symmetric and its density p(r) oc r" 1 55 from the 
shock to infinity, the normalization of the spectral fit gives an 
electron density of n e = 1 .06 ± 0.03 x 10" 3 cm" 3 at a radius of 
276 kpc (1 '.72). 

Using these parameters, the radius of the shock is 158 kpc, 
the time since the outburst is f s = 5.9 x 10 7 y and the total en- 
ergy of the outburst is £ s = 3 x 10 61 erg. This energy is simi- 
lar to the lobe enthalpy ( Gizani & Leahv 2004), as expected if 
the lobes drive the shock. The main source of uncertainty in 
the age of the outburst (~ 10%) is due to the uncertainty in the 
preshock temperature. The shortcomings of the model do give 
rise to systematic uncertainty in the s hock energy, but th is is 
unlikely to be more than a factor ~ 2 ( Nulse net al 120051) . 

In the temperature range 1.6 - 10 keV, the Chandra count 
rate in the band 0.6 - 7.5 keV is very insensitive to gas tem- 
perature, varying ±3.3% about its mean over the whole range 
of temperature, for a fixed emission measure. This enables us 
to deproject the gas density with reasonable accuracy, despite 
poor knowledge of the gas temperature (doubling the abun- 
dance to 1.0 would reduce the electron density by ~ 7%). 
A dep rojection was done, using the method of Nuls en et alJ 
(2005), with the gas temperature fixed at 4 keV and other pa- 
rameters as above. The resulting electron density profile is 
shown in Fig.0] together with the electron density profile ob- 
tained from the Mach 1.65 shock model. The results agree 
well with the model, clearly showing the density jump at the 
shock. The failure of the model for r < 25 kpc is a numeri- 
cal artifact, but other shortcomings are expected to make the 
model inaccurate at small radii. 

In the models, adiabatic expansion limits the size of the re- 
gion where the temperature of the shocked gas exceeds that of 
the unshocked gas. Nevertheless, the strength of the shock in 
Hercules A makes it a good candidate for detecting the tem- 
perature rise due to the shock. For the Mach 1.65 model, after 
projection onto the sky, the emission measure weighted tem- 
perature exceeds the preshock temperature by at least 20% for 
110 kpc < r < 150 kpc. With other fit parameters as above, 
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a spectrum extracted from this region (~ 900 0.6 - 7.5 keV 
source counts) gives a temperature of kT = 6. (90%), 
in reasonable agreement with the model. 



5. DISCUSSION 

The mean power of the outburst in Hercules A, P s = E s /t s = 
1.6 x 10 46 ergs _1 , is two orders of magnitude larger than 
the total power radiated from the region where the cooling 
time is shorter than 10 10 y. Hercules A joins a small col- 
lection of cooling flow clusters known to have large-scale 
sho cks driven by an outburst from an AGN at the cluster cen- 
ter JFabian et al.ll2003t iForman et alJl2005t FlVlcNama ra et alJ 
120051 iNulsen et alJ 1200511 Three of these systems, Hy- 
dra A, MS0735. 6+7421 and Hercules A, have outburst en- 
ergies of > 10 61 erg. The outburst in Hercules A cur- 
rently has the strongest shock and its total energy is the sec- 
ond largest known ( MS0735. 6+7421 is twice as energetic, 
McNamara et al.ll2005l) . Along with the other systems, it has 
important implications for the energetics of cooling flows, the 
preheating of clusters, the interaction of radio sources w ith the 
ICM and the growth of n uclear black holes dMcNamara et alJ 
120051 INulsen et all2005l) . 

If the outburst is powered by accretion onto a black hole, 
then the outburst energy is E s = eM s c 2 , where the mass M s 
was accreted to fuel the outburst and e is the efficiency of 
jet energy production by the black hole. Unless e > 10%, 
the mass swallowed by the black hole exceeds 1.7 x 10 8 M Q 
to fuel this outburst. If this mass was swallowed in a time 
comparable to the age of the outburst, t s ~ 6 x 10 7 y, the mass 
increase is hard to reconcile with a tight correlati on betwee n 
bulge properties and black hole mass (Gebhardt et al. 2000), 
unless the black hole is very massive indeed. 

In our Mach 1 .65 model, the shock inverts the entropy pro- 
file of gas inside 58 kpc (22"), creating a buoyant bubble at 
the cluster center that would then rise. A large bubble can rise 
at a significant fraction of the sound speed (IChurazov et al.l 



2001 ), but always more slowly than the shock front. Although 
the shock model is not expected to match reality closely, the 
cavities in Hercules A are comparable in size to the entropy 
inversion of the model, suggesting that this is how they were 
formed. This would explain the lack of radio emission from 
the cavities. 

6. CONCLUSIONS 

Analysis of a Chandra X-ray image of the Hercules A clus- 
ter shows that it has cavities and a shock front associated with 
the powerful radio source. Unusually, the cavities show no 
clear connection to the radio source. The shock front is elon- 
gated in the direction of the radio lobes and appears to be its 
cocoon shock. Fitting a simple hydrodynamic model to the 
surface brightness profile gives a Mach number for the shock 
of ~ 1.65. The age of the outburst that drove the shock is 
5.9 x 10 7 y and its total energy is 3 x 10 61 erg. The depro- 
jected density profile is consistent with the shock model and, 
in particular, with the density jump at the shock. Within the 
limits of the spectroscopic data, the temperature jump is also 
consistent with the shock model. 

The shock outburst is highly significant for the energetics 
of any cooling flow in Hercules A and for the cluster as a 
whole. The mean mechanical power of the outburst ~ 1.6 x 
10 45 erg s" 1 , well in the range of quasar luminosities. The 
black hole in the AGN that drove this outburst grew by, at 
least, 1.7 x 10 s M Q during the outburst. 
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tional Laboratory. 
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